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In the framework of four-band envelope-function formalism, developed earlier for spherical semi- 
conductor nanocrystals, we study the electronic structure and optical properties of quantum-confined 
lead-salt (PbSe and PbS) nanowires (NWs) with a strong coupling between the conduction and the 
valence bands. We derive spatial quantization equations, and calculate numerically energy levels 
of spatially quantized states of a transverse electron motion in the plane perpendicular to the NW 
axis, and electronic subbands developed due to a free longitudinal motion along the NW axis. Using 
explicit expressions for eigenfunctions of the electronic states, we also derive analytical expressions 
for matrix elements of optical transitions and study selection rules for interband absorption. 

Next we study a two-particle problem with a conventional long-range Coulomb interaction and 
an interparticle coupling via medium polarization. We derive analytical expressions for an effective 
direct Coulomb coupling and an effective coupling via medium polarization averaging corresponding 
coupling energies over the fast transverse motion of charge carriers, and then compute numerically 
the effective couplings for the lowest-energy electron-hole pair in a PbSe NW of the radius R — 5 ran 
in vacuum. The obtained results show that due to a large magnitude of the high-frequency dielectric 
permittivity of PbSe material, and hence, a high dielectric NW/vacuum contrast, the effective 
coupling via medium polarization significantly exceeds the effective direct Coulomb coupling at all 
interparticle separations along the NW axis. 

Furthermore, the strong coupling via medium polarization results in a bound state of the longi- 
tudinal motion of the lowest-energy electron-hole pair (a longitudinal exciton) , while fast transverse 
motions of charge carriers remain independent of each other. For a PbSe NW of the radius R = 5 nm, 
the binding energy of the longitudinal exciton is found to be about 77.9 meV that approximately two 
times smaller than the energy of spatial quantization of the lowest-energy electronic states. Thus, 
the strong interparticle coupling via medium polarization in quantum-confined lead-salt NWs signif- 
icantly modifies the single-particle electronic spectrum, and could result in essential modifications 
such Coulomb phenomena as impact ionization. Auger recombination and carrier multiplication. 

PACS numbers: 73.21. Hb, 73.22.Dj, 78.67.Lt 



I. INTRODUCTION 

Knowledge of the electronic states is needed for stud- 
ies of many physical properties such lower-dimensional 
quantum-confined semiconductor systems as nanocrys- 
tals (NCs) and nanowires (NWs). The electronic states 
of NCs have been studied analytically and numerically 
for both narrow-gap^ and wide-gap^ semiconductor ma- 
terials with different bulk band structures, while the 
structural, electronic, and optical properties of quantum- 
confined NWs have been computed mainly with ab ini- 
tio methodsj^i^i^i^ii such as density-functional theory, 
pseudopotentials^ and tight-binding methodsj ^°i^^'^^i^^ 
Theoretical studies of optical properties of NWs made 
of wide-gap materials, such as GaAs, CdSe, and InP, 
have been carried out^ in the framework of "particle-in- 
a-box" model not accounting for an interband coupling 
between the conduction and the valence bands. 

In this paper, in the framework of four-band envelope- 
function formalism, we first obtain analytical expressions 
for eigenfunctions of the electronic states of quantum- 
confined narrow-gap lead-salt (PbSe and PbS) NWs with 
a strong coupling between the conduction and the valence 
bands, and then derive the electronic structure imposing 
the boundary condition of vanishing envelope functions 
on the NW interface. Thus, we generalize an approach 
developed by Kang and Wise^ for spherical lead-salt NCs 



to the case of NWs of the cylindrical geometry. However, 
such a generalization is not trivial, because it requires 
new mathematical constructions, specific for the cylin- 
drical geometry, while in the case of the spherical geom- 
etry, conventional mathematical technique developed in 
the Dirac theory of relativistic electronic is employed. 

Making use of explicit expressions for electronic eigen- 
functions, we derive analytical expressions for matrix el- 
ements of the operator e • p (where e and p are the light 
polarization and the electron translation-momentum vec- 
tors, respectively), which determine optical absorption, 
and study selection rules for interband absorption. 

Then, we study a long-range Coulomb interaction in 
NWs. As in the case of NCsr^iii it contains a conven- 
tional Coulomb interaction, a coupling between a charge 
particle and a medium polarization created by this parti- 
cle itself, and finally an interparticle coupling via medium 
polarization, i.e. a coupling between a particle and a 
medium polarization created by the other particle. 

The interaction of a charge particle with its own "im- 
age" results in a particle repulsion from the NW interface 
that effectively reduces the NW radius, and thus slightly 
modifies the single-particle electronic spectrum. 

Since a transverse motion of charge carriers in the 
plane perpendicular to the NW axis in sufficiently long 
NWs of the length L > i?, where R is the NW radius, is 
obviously much faster than a longitudinal motion along 
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the NW axis, we derive analytical expressions for an ef- 
fective direct Coulomb coupling and an effective coupling 
via medium polarization averaging the energy of electron- 
hole (e-h) coupling with the eigenfunctions of the trans- 
verse motion of charge carriers. 

Further numerical calculations for a PbSe NW of the 
radius R = b nm in vacuum show that owing to a large 
magnitude of the high-frequency dielectric permittivity of 
PbSe material [Koo(PbSe) = 23], and hence, high dielec- 
tric NW/vacuum contrast, the effective direct Coulomb 
interparticle coupling is much weaker then the effective 
interparticle coupling via medium polarization at all in- 
terparticle separations along the NW axis. 

Solving numerically an effective eigenvalue problem for 
a relative longitudinal motion of the lowest-energy e-h 
pair, we found a bound state of the longitudinal motion 
(a longitudinal exciton) , and compute its binding energy, 
which is found to be about 77.9 meV that is approxi- 
mately two times smaller than the energy of spatial quan- 
tization of the lowest-energy electronic states. Therefore, 
in sharp contrast to the case of NCs, the Coulomb inter- 
action results in essential corrections to the spectroscopy 
of quantum-confined lead-salt NWs, and could signifi- 
cantly modify such Coulomb phenomena as impact ion- 
ization. Auger recombination and carrier multiplication. 



II. FOUR-BAND ENVELOP FUNCTION 
FORMALISM FOR CYLINDRICAL GEOMETRY 

In the four-band envelope-function formalism, total 
electronic wave functions in lead-salt semiconductor ma- 
terials are written as a product of the four (for two possi- 
ble directions of the electron spin in the conduction and 
the valence bands) band-edge Bloch functions \ui) in the 
conduction {i — 1,2) and the valence {i = 3,4) bands, 
and four-component envelope functions J-i, 



i=4 



(2.1) 



Boundary conditions in quantum-confined NCs and 
NWs, are imposed on envelope functions, which are found 
as solutions of the eigenvalue problem 



HT = ET 



(2.2) 



with the eigenenergy E and the Hamiltonian of bulk lead- 
salt materials in a spherical approximation^ 
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Here p ~ — iV is the wave vector operator applying to 
the envelop functions, cr = (ax, cry, a z), 
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are the Pauli matrices, rj = ■:^P (where mg is the free 
electron mass, and P is the Kane momentum) is a pa- 
rameter of the interband coupling. 
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are the operators of bare (i.e. in the absence of the in- 
terband coupling) electron energies in the conduction (c) 
and the valence (u) bands, E^ is the energy gap, and ra^.v 
are the effective electron masses in the bands. 

In the cylindrical coordinates {r,(j),z) = {r,z), where 
the Z axis is directed along the NW axis, and r is the 
radius-vector in the XY plane perpendicular to the Z 
axis, it is convenient to separate the transverse and lon- 
gitudinal motions and to rewrite the Hamiltonian as 



where 



and 



Et — El^y + Hz , 
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Here, the total wave vector operator p is represented as a 
sum of the wave vector operator q = — of the trans- 
verse motion and the wave vector operator 
the longitudinal motion. 

It is convenient to find first solutions of an auxiliary 
eigenvalue problem with the Hamiltonian of the trans- 
verse motion H^y, and to study then the eigenvalue prob- 
lem (2.2) with the total Hamiltonian H . 



III. AUXILIARY EIGENVALUE PROBLEM 



The auxiliary eigenvalue problem 



(3.1) 



for a four-component bispinor iIj = ^ ^ J ' where and 
X are two-component spinors, takes the form 

[eo - f-ciq)] If ^ -nia ■ q)x, (3.2a) 
[eo + ec(q)] X ^ vifT ■ qW- (3-2b) 

In the absence of the interband coupling (particle-in- 
a-box model), 77 = 0, solutions of Eqs. (3.2) are easily 
found to be 



where J„i, {qr) are the Bessel functions, 
1 A . / 



(3.3a) 
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are the spinors corresponding to two possible projections 
of the electron spin on the Z axis, and to; = 0, ±1, . . . are 
eigenvalues of the orbital angular momentum operator of 
the transverse motion l^rxq = —i(^rx ^) directed 
along the Z axis. 

Taking into account also a free longitudinal motion 



with the eigenfunction 
particle- in-a-box model 



the eigenenergies in the 
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correspond to states in the conduction (+) and the va- 
lence band (— ) with the total wave vector p = y/q^ + k^, 
where q and k are wave vectors of the transverse and the 
longitudinal motions, respectively. 

The wave functions (3.3a) describe states with given 
z-projections of both spin and orbital angular momen- 
tum of the transverse motion. However, because of 
the interband spin-orbit coupling term in the Hamilto- 
nian H^y, these projections are not conserved separately, 
and only the projection of the total angular momentum 
rrij = mi -\- Sz is conserved. 

From functions e™'*^ and spinors £,a (ct = ±) one can 
construct two polar angular spinors: 



^ruj {4>) 
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(3.4a) 



where = r /r = cos (j)- sin <p ■ e 
unit vectors along the X and Y axis, and hence 



and By are the 
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e"!' 



The relations (3.6) play an important role in further 
computations, because they allow one to split the sys- 
tem of equations (3.2) into independent equations for the 
bispinor components Lp and x- 

As in the case of the spherical geometry;^ one can 
now construct two different bispinors with a given total 
angular momentum projection nij and uncertain values 
of the orbital angular momentum and spin projections: 



^/m,(r)an,(?!') 



and 



^ _ /C/™;(r)e„^(0) 



(3.7a) 



(3.7b) 



where the coefficients and radial functions /,„, and /„j 
should be found from solutions of the eigenvalue problem. 

Inserting the expressions (3.7) into Eqs. (3.2) we find 
after tedious but straightforward computations: 
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and 
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rrii = mj + 



(3.4b) 



which are analogous, in some sense, to the angular 
spinors^^ in the spherical geometry. 

The polar angular spinors are eigenstates of all three 
angular operators: the spin operator — with 
eigenvalues Sz = the total angular momentum op- 
erator j = I + Sz with eigenvalues rrij, and the oper- 
ator I with two possible eigenvalues m/ = rrij — i and 
m'l = rrij -|- ^ at a given value of the total angular mo- 



mentum projection. Moreover, they are orthonormal. 
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m-jm ■ -i 



m-jm ■ -i 



(3.5a) 
(3.5b) 



and orthogonal to each other, 8,^ i^m- = 0, owing to 

j 

orthogonality of the spinors ^+ and ^_ . 

The introduced polar spinors are related to each other 
by the expressions 



^raj — z(^T ■ e^)0^^. , 



(3.6a) 
(3.6b) 



and 
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~ Cg/ m, ir)i^m, (0) 
-\/eO + eg/m;(»')0mj (</>) 



\/eo + eqf ni[ {r)<dmj (</>) 

V^O - £qfni, {r)^nij (0) 



\/ £o - eg/ r«; {r)Qm, (0) 



(3.8b) 



(3.9a) 



(3.9b) 



in the conduction (+) and the valence (— ) bands, which 
describe states with the eigenenergy 



£o = 



jy2g2 



(3.10) 



in the conduction band, and — £o in the valence band. 
To simplify the above expressions, we used a "mirror" 
symmetry of the conduction and valence bands in PbSe 
and PbS materials, and set rric = rriy = m. Therefore, 

here and hereafter eg = ^ + -2^- 
The radial functions 



fm, {r) 
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where (Ar) are the modified Bessel functions, A 
y?TAj, and Aq 



2m 



+ rf, vanish on the NW 
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interface (r = R) at any wave vector q. The coefficients 
Nfni are found from the normalization condition 



/' 



With normahzed /-functions, the bispinors are obviously 
normalized to unity. 

The condition of self-consistency of simultaneous van- 
ishing both spinor components of the bispinors on the 
NW interface results in the spatial quantization equa- 
tions for the transverse wave vector: 



^£0 -I- eg Jmi {qR) _ ygp - eg Jm[{qR) 

Vex - £0 Imi (Ai?) -v/cATeo Im[ (Ai?) 

for the bispinors and and 

V£o + eg Jm[ {QR) _ ygp - eg Jm, [QR] 
\/eA - £o Im[{>^R) \/eA + £o Imi{>^R) 



(3.12a) 



(3.12b) 
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for the bispinors ^_ and Here, eA = 

EA = £q + ^^'^ '^'^ ^''C two notations to emphasize 

the existence of two distinct sets of spatially quantized 
wave vectors, g„ and Qn (n = 1, 2, . . .), corresponding to 
the same total angular momentum projection mj. 

As it must be expected, ^t-bispinors are orthogonal to 
^-bispinors due to orthogonality of the fl and 6 polar 
angular spinors, while orthogonality of bispinors corre- 
sponding to states in different bands, and $-|-, 
with the same mj is provided by the structure of 
energy factors in the expressions (3.8) and (3.9). 

The bispinors describing electronic states in the con- 
duction band contain a contribution of the valence band, 
and vice versa. Correspondingly, the total wave functions 
(2.1) contain contribution of the band-edge Bloch func- 
tions of the conduction and the valence bands. While, in 
the absence of the interband coupling, the second compo- 
nents corresponding to the valence band in ^+ and $+ 
bispinors and the first components corresponding to the 
conduction band in and bispinors vanish due to 
vanishing the factor ^£o — at ?7 — > 0. 



IV. LONGITUDINAL MOTION 

The longitudinal motion along the Z axis additionally 
mixes quantum states in the conduction and the valence 
bands. Solutions of the eigenvalue problem with the total 
Hamiltonian 



(4.1) 



can be found as a linear superpositions of the eigenstates 
of the transverse motion Hamiltonian H^y 

Since the interband coupling term in the Hamiltonian 
Hz mixes quantum states in the conduction and the va- 
lence bands describing by different bispinors \E' and 



one can look for solutions of the eigenvalue problem (4.1) 
in the form 

F = (y4*+ -^S$_)e^'=^ (4.2a) 
G = (C*_ -h ^>$+)e''=^ (4.2b) 

where the coefficients should be found from solution of 
the eigenvalue problem. Inserting these expressions into 
Eq. (4.1) we finally find: 



V = ( \/E + £*+ - VE - £$_") 

V2E y / 



(4.3a) 



;+ = (yE^^_ + e'''' (4.3b) 

for states in the conduction band with the eigenenergy 
E = +^el + rpp\ (4.4) 

and 

F_ = (^Ve^^+ + /BTi$-) e'''' (4.5a) 



Ve + s-^- - Ve~- 



£$+) 



(4.5b) 



for states in the valence band with the eigenenergy 
Here £ = 



/ei + : 



> m% ^P = ^ + ^> and = 
It is easy to see that the longitudinal motion does not 
modify the spatial quantization equations (3.12), because 
it mixes bispinors with the same sets of the transverse 
wave vectors g„ and Qn- 

Thus, eigenstates of the total Hamiltonian H in the 
cylindrical geometry are characterized by the projection 
of the total angular momentum m^- = ±i, ±|, . . . on the 
Z axis and the continuous wave vector of the longitudinal 
motion fc, while the spatially quantized wave vectors qn 
and Qn of the transverse motion are determined by Eq. 
(3.12a) for the bispinors F±, and by Eq. (3.12b) for the 
bispinors G± . The bispinors are orthogonal to each other 
and normalized to 2TTS{k — k'). 

Since there are two kinds of bispinors {F and G), the 
number of quantum numbers, which characterize eigen- 
states of the Hamiltonian H, is equal to that in the 
absence of the interband coupling (the particle-in-a-box 
model), as it must be expected. 

Finally, corrections related to mixing of quantum 
states in the conduction and the valence bands due to 
the longitudinal motion are essential only at sufficiently 
large vectors k comparable to magnitudes of the trans- 
verse motion wave vectors q and Q. While the densities 
of electronic states per unit length of a NW exhibit ob- 
viously the Van Hove singularities at fc ^ 0: 
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where is the number of states, and £0(9) are the 
subband-edge energies determined in Eq. (3.10). 
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FIG. 1: Band-edge energy levels {k = 0) in the conduction 
band in the absence of the interband coupling. 
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V. ELECTRONIC STRUCTURE 

For further comparisons, we consider first the elec- 
tronic structure in the framework of particle-in-a-box 
model, in which the boundary condition ipir — R) — 
for the wave functions (3.3a) results in a simple spatial 
quantization equation 



J™,(gi?) = 0. 



(5.1) 



The first three zeroes jm,,n for m; = 0, ±1, ±2, where the 
second index n = 1,2,... numerates zeroes of the Bessel 
functions are approximately found to be: 

ja,i = 2.4 jo, 2 = 5.5 jo.3 = 8.7 
jis = 3.8 ji,2 - 7.0 ii^3 = 10.2 

32,1 = 5.1 j2,2 = 8.4 j2,3 = 11-6 

Then, the subband-edge energies (fc — 0) corresponding 
to states with the orbital angular momentum projection 
±mz are given by: 



^m.i .71. ^ 
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2 .-2 

mi ,n 



2mi?2 



(5.2) 



The subband-edge states are additionally degenerate 
with respect to two possible projection of the electron 
spin on the Z axis. Thus, we find for the energies of 
subbands with the orbital angular momentum projection 
±771; and the continuous wave vector k 



2m 



(5.3) 



where the signs -I- and — correspond to the conduction 
and the valence band, respectively. The subband-edge 
structure in the conduction band is shown in Fig. 1. 

The interband coupling completely lifts the degener- 
ation of subband-edge states. The results of numeri- 
cal calculations of the electronic structure in the con- 
duction band for a PbSe NW of the radius R = 5 nm 
are presented in Table I and illustrated in Fig. 2. For 



FIG. 2: Spatially quantized electronic states described by the 
bispinors ^'+(g) and $+(Q) in the conduction band, and mir- 
ror symmetric states in the valence band described by the 
bispinors "l>-(g) and $_(Q). 



State (mi) 





+ 1 -1 +2 


-2 


Energy (eo) 


0.22 


0.32 0.33 0.43 


0.45 


State (m'l) 





+ 1 -1 +2 


-2 


Energy (eo) 


0.25 


0.37 0.37 0.49 


0.49 



TABLE I: Energies of the spatially quantized states de- 
scribed by bispinors 'i'+(q) (mi = 0, ±1, ±2) and <I?+(Q) 
{m'l = 0, ±1, ±2) in the conduction band. The states mj — ±1 
and m'l = ±2 are actually not degenerate, but the energy dif- 
ferences are of the order of a few meV only. 



numerical calculations we adopted the following param- 
eters of PbSe material from Ref. [1]: Eg = 0.28 eV, 
rric = rriy = 0.20mo, and 7] = 0.31 eV-nm~^. 

Although the orbital angular momentum projection is 
not a good quantum number, it is still convenient to char- 
acterize quantum states by the orbital angular momen- 
tum projection of the first (second) spinor component of 
bispinors in the conduction (valence) band. The projec- 
tion of the total angular momentum rrij is determined by 
the expressions rrij = mi + ^ and mj = m'l — ^. Thus, 
the lowest (upper) state in the conduction (valence) band 
corresponds to mj — ^. 

Finally, Table II summarizes magnitudes of the effec- 
tive energy gap in PbSe NWs of different radii computed 
in the framework of four-band envelope function formal- 
ism (the second row) and in the particle-in-a-box model, 
77 = 0, (the third row). It is easy to see that the interband 
coupling results in significant size-dependent corrections 
to the electronic structure computed within the frame- 
work of the particle-in-a-box model. Note that the inter- 
band coupling decreases the wave vector of the lowest- 
energy electronic states, but increases the effective en- 
ergy gap, while in lead-salt NCs the interband coupling 
reduces^ E* 
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R (nm) 4 5 6 7 8 

El (cV) 0.51 0.44 0.40 0.37 0.35 
K (cV) iji = 0) 0.42 0.37 0.34 0.32 0.31 



TABLE II: The effective energy gap of PbSe NWs of different 
radii calculated in the framework of four-band envelope func- 
tion formalism (the second row) and particle-in-a-box model, 
77 = 0, (the third row). 

VI. INTERBAND OPTICAL TRANSITIONS 

The strength and selection rules of optical transitions 
in NWs are determined by the matrix element 

M = {iP\e-p\iP), (6.1) 

where e is the photon polarization vector and are the 
total electronic wave function defined in Eq. (2.1). 

For the case of intorband transitions, the matrix ele- 
ment Mcv within a given L- valley of lead-salt material is 
derived by the method described in Ref. [1] as 

Mcv = {e-z)Pi j drdzJ^l{r,z){ax®cr^)J^y{r,z) 

+ j drdzTl{r,z){e-p)Ty{r,z), (6.2) 

where ® stands for the direct product. To simplify fur- 
ther expressions, we rewrite Eq. (6.2) as 

Mcv = (e • z)Pi{TcWx ® (^z\J^v) + {J^c\e- ■ p\Tv) 

In the first term, z stands for one of four equivalent (111) 
directions in the face-centered cubic lattice, and Pi is the 
matrix element of tlw^ longitudinal Kane momentum be- 
tween the band-edge Bloch functions. This term becomes 
isotropic as result of summing over all four equivalent L 
valleys. 

Inserting into the first term of Eq. (6.2) the expressions 

for and bispinors, and using results of computa- 
tions of matrix elements in Appendix A, we find 

M^fIf. = (e--2)-Pz(^V,m^,fc'|o-^®0"^|-P-,mi,fc) 

= -2w{e-z)Pi^S{k'-k)5m'.m,, (6.3) 
where the ratio e/E can be written as 

This expression is valid only for transitions with the same 
vectors of the transverse motion q'^ = Qn corresponding 
to m!j = TTij. Transitions with ^ qn are also allowed, 
but their matrix elements are smaller due to reducing 
the radial functions overlap at different qn- The ma- 
trix elements for the "direct" transitions from G- to 
subbands differ in sign only, Mq^ q_ = —Mp^p_, while 



the matrix elements M^^q^ and M^^p^ for allowed 
{m'j = rrij) "indirect" transitions are small at fc <C Qn, 
because they are proportional to a small factor r]k/E*. 

It should be emphasized that S{k' — k) function, which 
expresses conservation of the Z-component of the wave 
vector in the system electron plus photon, appears in 
Eq. (6.3a) because we neglect the small photon wave 
vector. If the photon wave vector is taking into account, 
(5-function takes the form 5{k' — k — Kz), where Kz is the 
projection of the photon wave vector on the Z axis. 

Thus, the first term of the matrix element is mainly 
determined by the longitudinal Kane momentum of bulk 
material, while envelope functions determine the ratio 
e/E (for transitions with q'^ = qn), which equals unity 
at the subband edges, and is reduced when k grows up. 
Taking into account Van Hove singularities in the density 
of electronic states at A; —> 0, we should conclude that 
the subband-edge absorption significantly exceeds that 
at finite k. 

On the contrary, the polarization-dependent second 
term of the matrix element in Eq. (6.2) (which is obvi- 
ously absent in the particle-in-a-box model) is completely 
determined by envelope functions. The matrix elements 
for the direct transitions Mp^p_ = (-F+le • p|F_) and 

M^^^Q_(G+|e.p|G_) are proportional to ryfe/E"*. There- 
fore, their contribution to the transition strength is small 
in comparison with that of the matrix elements Mp^p_ 

and Mq^q_. While matrix elements for indirect transi- 
tions are found to be at small k: 

mP^G. = (*+,m^|e-q|*-,™,), (6.4a) 
^gIf. = (*+,m^.|e.q|$_,™,). (6.4b) 

They do not vanish only for transitions with m'j = mj±l. 

Finally, matrix elements of the operator Sz • kz van- 
ish owing to orthogonality of states in the valence and 
the conduction bands, and absorption of light polarized 
along the NW axis is determined by the matrix elements 
and M^Iq_ only. 

VII. COULOMB INTERACTION AND 
LONGITUDINAL EXCITONS 

Now we are able to consider a two-particle problem 
in quantum-confined lead-salt NWs. To solve the two- 
particle problem is convenient to write the NW Hamil- 
tonian in the absence of an interparticle coupling in the 
second-quantization representation as 

^0 = E / Hik)cgik) + hl{k)h,{k)] . (7.1) 

Here, the operators cJ(A;) [cg(fc)] and /iJ(A;) [hq{k)] create 
(annihilate) an electron and a hole with a wave vector 
in a subband with a spatially quantized wave vector 
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q in the conduction and the valence band, respectively, 
and the summation is carried out over all the subbands. 
In what follows, we restrict our consideration to the case 
when both electron and hole belong to the lowest-energy 
subbands F, i and F_ i . That allows one to omit the 

^'2 '2 

summation over other subbands in Eq. (7.1). 

For a sufficiently slow longitudinal motion of charge 
carriers with fc <C <?, where now q is the wave vector of 
the transverse motion in the lowest-energy subbands, the 
expression for the eigenenergy (4.4) is expended as 



Eq{k) ~ eoiq) + — = -E. 



2 1,2 



where 



2m,(g) 2 s 2m,{q)' 



1 E* 

mM ^ 



(7.2) 



is the effective mass of the longitudinal motion. Intro- 
ducing the Fourier transformed charge carrier operators 



ikz 



one can rewrite the free Hamiltonian (7.1) as 



2 ^2 



dzhh\{zii 



Cq{Ze) 



2mz dzf^ ' 



hg{zh) (7.3) 



At fc ^ q, the energy factors \JE — e in Eqs. (4.3) 
and (4.5) are also small, and the terms proportional to 
\/E — e can be omitted. Then, the total bispinors F± 
differ from and <I>-bispinors in the wave function of 
a free longitudinal motion only, i.e. _F1|_ ~ 4'+ exp (ikz) 
and -F- ~ exp (i/cz). Correspondingly, the operator 
of e-h coupling is written as 

^ = J dZedZhc\{Ze)Cq{Ze)V{Ze- Zh)h\{zh)hq{zh), 

(7.4a) 

where the effective coupling energy is given by 

y = J dr-edr/i^'^ i(re)^'+ i(re)C/(re,2:e;r/i,2,,) 



x$^_ i(r,,)$_ 1 (r/i). 



(7.4b) 



with the coupling energy U determined in Eqs. (B5). 

Then, the eigenvalue problem for a relative motion of 
the lowest-energy e-h pair with the total Hamiltonian 
H = Ho + V reads 



rriz d(^ 



— :j7^MC) + [seh - v{C)] MO = 0, (7.5) 



where C = Ze-Zh, e^h = - El , and E^^ is the energy 
of the pair. Due to the strong interband coupling, the 
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FIG. 3: The self-interaction energy vs. normalized radial 
coordinate for a PbSe NW of _R = 5 nm; the constant term 
f4eif(f' = 0) = 57.47 meV is extracted. 



effective mass of the longitudinal motion, niz, in a PbSe 
NW of the radius of 5 nm, is found to be about a half of 
the effective electron mass m, ~ 0.51m. 

A charge carrier confined in a cylindrical NW of the 
dielectric permittivity placed in a host of the per- 
mittivity Kh creates a medium polarization (an "image" 
charge) that results in an interaction between the charge 
and the image charge. The expression for their interac- 
tion energy, Usqu, is derived in Appendix B [Eq. (B4)], 
and the results of numerical calculations for a PbSe NW 
of the radius R = 5 nm in vacuum are plotted in Fig. 
3. This single-particle term must be included in the 
Hamiltonian (2.3) determining the single-particle elec- 
tronic spectrum. Since Uscuir) grows up with r and di- 
verges at r = i?, the self-interaction repulses a charge 
carrier from the interface that obviously results in effec- 
tive decreasing the NW radius and increasing (decreas- 
ing) electronic energies in the conduction (valence) band. 
However, the magnitude of Usc\f{r) reaches magnitudes 
of confinement energies only at r ~ 0.99i?. Therefore, 
corrections to the single-particle spectrum are inessential 
and can be omitted. 

As in the case of NCsjA^'l^ the energy of total e-h in- 
teraction U (B5a) is separated into the energy of direct 
Coulomb coupling, Uc [see Eq. (B5b)], and a term cor- 
responding to the interaction energy between one charge 
carrier and a medium polarization created by the second 
one, Upo\ (B5c). For the lowest-energy e-h pair, the direct 
Coulomb (Vc — {Uc)) and polarization (V^oi = (f^poi)) 
parts of the total effective e-h coupling are plotted in Fig. 
4 as functions of the modulus of e-h separation along the 
NW axis normalized to the NW radius. It is easy to see 
that due to a large magnitude of the dielectric permittiv- 
ity of PbSe material [^(PbSe) — 23], and hence, high di- 
electric NW/vacuum contrast, the effective coupling via 
medium polarization Vpoi essentially exceeds the effective 
direct Coulomb coupling Vc at all e-h separations along 
the NW axis. Therefore, the latter can be omitted in the 
eigenvalue problem (7.5). 

To avoid too tedious computations with the very com- 
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FIG. 4: Energies of the effective direct Coulomb electron- 
hole coupling Vc (top) and the effective coupling via medium 
polarization Vpoi (bottom) vs. the normalized modulus of e-h 
separation ~ \ze — Zh\/R along the NW axis for a PbSe 

NW of i? = 5 nm. 

plicated function Vpoi(C) derived in Appendix C and plot- 
ted in Fig. 4 (bottom), it is convenient to replace V^oi(C) 
by the function 

V{C/R) = , " (meV), (7.6) 

V(C/i?)2 + 62 

which, at a « 392.7 meV and b k, 3.3, well interpo- 
lates the function Ipoi(C)- Then, numerical computa- 
tions show that the lowest-energy electron and hole form 
a bound state (a longitudinal exciton) with the exciton 
binding energy about 77.9 meV, Seh — —77.9 meV, while 
fast transverse motions of charge carriers remain inde- 
pendent of each other. The excitation energy of the lon- 
gitudinal exciton, Sg/i = £"* + Seu , is then estimated to 
be about 362 meV, and hence, it is placed approximately 
in the middle point between the NW effective energy gap 
i?* = 440 meV and the energy gap — 280 meV of 
bulk PbSe material. 

Finally, the size of the longitudinal exciton is estimated 
to be about lOi? — 50 nm that justifies the approxima- 
tion of a slow longitudinal motion in comparison with 
the transverse motion of charge carriers. The longitu- 
dinal exciton size is comparable to the Bohr radius of 
Wannier-Mott exciton, as = 46 nm, in bulk PbSe ma- 
terial. However, such large longitudinal excitons can be 
destroyed by charge carrier scattering on impurities and 
other imperfections in NWs.-S Therefore, the problem of 



longitudinal excitons and their role in the spectroscopy 
of quantum-confined lead-salt NWs requires further the- 
oretical and experimental studies. 

VIII. CONCLUSIONS 

In conclusion, in the framework of four-band envelope- 
function formalism we have studied the electronic struc- 
ture and optical properties of quantum-confined lead-salt 
NWs with a strong coupling between the conduction and 
the valence bands. Numerical calculations show that the 
interband coupling completely lifts the degeneration of 
electronic states and results in significant size-dependent 
corrections to the electronic structure computed in the 
framework of particle-in-a-box model. We have also de- 
rived analytical expressions for the matrix elements of 
the operator e • p, which determine optical absorption 
in NWs, and have studied selection rules for interband 
absorption. 

Finally, we have studied a two-particle problem with 
an effective long-range Coulomb coupling averaging cor- 
responding interaction energy over a fast transverse mo- 
tion of charge carriers. Numerical calculations show that 
due to a large magnitude of the dielectric permittivity 
of lead-salt materials, the effective interparticle coupling 
via medium polarization significantly exceeds the effec- 
tive direct Coulomb coupling. 

Furthermore, the strong coupling via medium polar- 
ization results in a bound state of the longitudinal mo- 
tion of the lowest-energy e-h pair (a longitudinal exciton) 
with the large binding energy comparable to the energy 
of spatial quantization of electronic states, and the size 
about 50 nm, which is comparable to the Bohr radius 
of Wannier-Mott exciton in bulk PbSe material. Thus, 
the strong Coulomb coupling in lead-salt NWs results 
in significant two- and many-particle corrections to the 
single-particle electronic spectrum, and could essentially 
modify such Coulomb phenomena as impact ionization. 
Auger recombination and carrier multiplication. 
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APPENDIX A: MATRIX ELEMENTS 

Here, we compute the matrix elements of the opera- 
tors (Jx 'S) Oz and e • p on the and <i>-bispinors, which 
are needed to compute the first and the second terms in 
the expression (6.2), and to derive the selection rules for 
optical transitions. 

Note that owing to orthogonality of the fi- and 0- 
spinors and the structure of the operator <Jx®<Jzi "diag- 
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onal" matrix elements vanish, i.e. 

^'^K«)(T^)* = $Ucr:r«)cr^)* = (Al) 

at any quantum numbers of the bispinors in both the 
conduction and the valence band. While "non-diagonal" 
matrix elements {'i>+\ax <E)(Jz\^-) and ((f>_|_|(Ta; (8) ctz I'l'-) 
are found to be: 



{'^+.mr\ax (g) crz|$_,™J 
($ + ,m'.|cr^«)Crz|^'-,m,) 



(A2) 

(A3) 



where we used Eqs. (3.5). Therefore, the matrix elements 
of the operator (g) tr^ are diagonal in the quantum num- 
ber rrij and correspond obviously to the direct interband 
transitions. Note that these simple expressions are de- 
rived for transitions between states with the same wave 
vectors = qn- 

On the contrary, owing to orthogonality of polar an- 
gular spinors Q and Q, non-diagonal matrix elements 
|e • p|$_^„j^.) and |e • pl^-.,™^) vanish at 

any quantum numbers. 



4'^(e • p)^ = $^(e • p)^' = 0. 



(A4) 



While diagonal matrix elements ('I'+.m' |e •p|^'__„i^) and 
i^+.m'. |e ■ p\^-^mj) do not vanish, and contribute to the 
strength of optical transitions. 



APPENDIX B: COULOMB INTERACTION 

To our best knowledge, an expression for the Coulomb 
interaction in NWs was derived for the first time in Ref. 
[14]. However, the derived expression contains some 
mathematical incorrectness, so that electrical potentials 
are complex rather than real as it must be. Therefore, 
here we present a correct solution of the Coulomb prob- 
lem, and derive an asymptotic expression for particle- 
particle coupling at interparticle separations along the Z 
axis, essentially exceeding the NW radius, \z ~ z'\ ^ B.. 

An electrical potential P{r,z) at the point with the 
coordinates (r, z) created by an electrical charge e posi- 
tioned at the point {r',z') inside a NW of permittivity 
Knw placed in a host medium of permittivity is deter- 
mined by the Poisson equation 



,AP(r, z) = -AneSir - r')S{z - z ) 



(Bla) 



inside the NW {r < R), and by the Laplace equation 

AP(r,z) = (Bib) 

outside the NW (r > i?), where A is the Laplace op- 
erator. Solving these equations with the conventional 
boundary conditions^ 



P{r = R-0) = P{r = R + 0), 

Kr,y,-^Pir ^ R-0) = Kh-^P(r = i?-f 0),(B2b) 
ar ar 



(B2a) 



we find for the potential inside the NW (r < R) 

Pin{r, z) = Pc{r, z) + Ppoi(r, z), (B3a) 

where 

p„(r z) = — ^ 

Hn^y/\r-r'\'^ + {z-z'f 

= ^ e™(*-*') r^cos[fc(z-z')] 

m— — oo ^ 

^ f I^{kr)K^{kr'), r < r' 
\ K^{kr)I^{kr'), r>r' 

is the Coulomb potential, and 

4e 



(B3b) 



Jm{<t>-(t>') 



m— — oo 



X 



f°° dk 

I — COs[/c(z-z')]Qm(/ci?) 

xlraikr)lraikr') (B3c) 



is the potential created by the medium polarization, or, 
in other words, the potential of an "image" charge. Here 



,ikR) = - 



K'^{kR)K„,{kR) 



K„,{kR)r„^{kR) - ^K'^{kR)I„^{kRy 



Im and Km are the modified Bessel functions, and 
rm{z) = fjrn{z),KUz) = ±Km{z). 

Note that integration over k in the Fourier integrals in 
Eqs. (B3), and hence, in all further expressions of this 
kind, is performed over positive values of k only, that 
eliminates imaginary values related to the imaginary part 
of the functions Km at negative argument. 

As in the case of nanocrystals )^^'^^ the interaction en- 
ergy of two particles with the charges e and e' in a NW 
is separated into a self-interaction energy, i.e. the energy 
of interaction of a charge with its own image: 



C/,eif(r) = iePpoi(r ^ r\ z ^ z') = — (l ^ ^ 



m— — oo 



dk 



J2 I ^QrnikR)llikr), (B4) 



and the energy of interparticle coupling 

U = Uc + C/poi. (B5a) 
Here, the first term is the direct Coulomb coupling energy 



Knw y/\r - r'p + (z - z')^ 
^ J- ^«n(4>-4>') r ^ COS [k{z-z')] 

Knw Jo 27r 

m— — oo ^ 



Imikr)Kmikr'), r < r' 
Km{kr)Im(kr'), r > r' 



(B5b) 
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while the second one is the energy of interparticle cou- 
pUng via a medium polarization, i.e. an interaction of 
one charge with the medium polarization created by the 
second one 



4ee' ( , Kh \ „,:m(0-0') /"^ * 



27r 



X cos \k{z - z')\Qra{kR)l^{kr)Im(kT'). (B5c) 

At large interparticle separations \z — z'\ 3> -R, the 
leading asymptotic term of the polarization part of the 
interparticle coupling is found to be 



ee' 1 



pol 



(B6) 



Kcff \Z - Z'\ ' 

where the effective permittivity is given by 
1 _ 1 1 

^cff '*h '^nw 

While the asymptotic of the direct Coulomb coupling 
reads 



Ur 



ee' 1 



(B7) 



APPENDIX C: EFFECTIVE ELECTRON-HOLE 
COUPLING 

Inserting into Eq. (7.4b) the expression for the energy 
of direct electron- hole Coulomb coupling (B5b) , we find 



for the effective Coulomb coupling 



4e2 r°° dk 

Vc(t) = / — cos(fcr) / xiX2dxidx2 

l^nw^ Jo Jo 



xD{xi)KQ{kxi)Io{kx2)D{x2), (Cla) 



where 



Di^) - ^ [(^0 + e,)f^,{x) + (£0 - e,)f,ix)] , 

and we introduced dimensionless variables xi = r^/R, 
X2 — rh/R, and t — \ze — Zh\/ R. Note that the angle in- 
tegration in Eq. (7.4b) selects the only nonzero term with 
m — from the sum over m in Eq. (B5b). Integrating — 
over k in Eq. (Cla), we finally derive 

e2 

^c(t) = ^ / xidxiX2dx2D{xi)D{x2) 

Knw^ Jo 
F ( ^ i- 1 ■ i'A'A 



^Jx\+ x'I+t'^ 



(Clb) 



where 2F1 is the hypergeometric function. 

Analogous computations for the polarization part of 
the effective e-h coupling result in 



T/ I \ 4e2 / \ dk 

r^nwJi- \ l^nw /JO 
2/ 



xgo(fc)J'(fc), 

where J{k) = xdxD{x)lQ{kx) 



(C2) 
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